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Abstract 

We study the impact of service-time distributions on the distribution of the maximum 
queue length during a busy period for the M-^ /G /\ queue. The maximum queue length is 
an important random variable to understand when designing the buffer size for finite buffer 
(M/G/l/n) systems. We show the somewhat surprising result that for three variations of 
the preemptive LCFS discipline, the maximum queue length during a busy period is smaller 
when service times are more variable (in the convex sense). 
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1 Introduction 



An important design issue for telecommunication systems and other applications is determining 
the buffer size when buffers are finite. We can better understand the effect of a particular buffer 
size by understanding the distribution of the maximum queue length during a busy period 
in an infinite buffer system. We give a characterization of the busy-period maximum queue 
length, M, for the queue for three types of preemptive LCFS (last-come first-served) 

disciplines: (i) preempted services are resumed when service recommences (LCFS-p-resume) , 

(ii) preempted services must be restarted from scratch when service recommences and a new 
service time is chosen from the service-time distribution (LCFS-p-repeat with resampling), and 

(iii) preempted services must be restarted from scratch when service recommences but the total 
service requirement for a given customer is the same each time it restarts its service (LCFS-p- 
repeat-without-rcsampling). These characterizations of M for each of the queueing disciplines 
allow us to show the effect of service-time distributions on M, as stated in (i)-(iii) below. For 
a fixed service discipline, let M and M' be the maximum number of customers during a busy 
period in two /G/1 queues with respective generic service times S and S', and with the 
same arrival rate A, and the same batch-size distribution. We assume that the distributions of 
S and iS" are such that the queues are stable. In this paper we show that the following relations 
hold; the definitions of the various stochastic orders can be found in the next section. 

(i) Under the LCFS-p-resume discipline, if S' <lt S, then M' <st M. 

(ii) Under the LCFS-p-repeat (with resamphng) discipline, if E(e~^^') > E{e~'^^), then 

M' <st M. 

(iii) Under the LCFS-p-repeat-without-resampling discipline, if S' <icv S, then M' <st M. 



A consequence of our results is the somewhat surprising conclusion that M will be stochas- 
tically smaller when service times are more variable (in the convex sense) under the preemptive 
LCFS disciplines. Miyazawa (1990) and Miyazawa and Shanthikumar (1991) show that for the 
finite-buffer M^/G/l/n queue under a non-preemptive discipline, the loss rate, i.e., the proba- 
bility that a random customer is lost, will be larger when service times are more variable in the 
convex sense. Our result relates to the loss rate, but the effect goes in the other direction. That 
is, we have that for preemptive LCFS disciplines, P(M > n) is smaller when service times are 
larger in the convex sense, where P{M > n) can be interpreted as the probability of at least 
one loss during a busy period in the /G/l/n queue. See also Chang, Chao, Pinedo, and 
Shanthikumar (1991). 
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For other results on the impact of the service time and batch size distributions on various 
performance measures of queueing systems, see, for example, Hordijk (2001), Makowski (1994), 
and Shanthikumar and Yao (1994), and the references therein. For other applications of the 
preempt-repeat service discipline, see, e.g., Adiri, Frostig, and Rinnooy Kan (1991), Birge, 
Frenk, Mittenthal, and Rinnooy Kan (1990), Cai, Sun, and Zhou (2004), and Cai, Wu, and 
Zhou (2004). 

The paper is organized as follows. We first recall some definitions of stochastic ordering in 
the next section. We then study M for each of the preemptive LCFS disciplines. Finally we 
provide some numerical illustrations of our results. 

2 Preliminaries 

Recall the following stochastic ordering relations for random variables X and Y . 

Definition 2.1 X is larger than Y in the stochastic sense, X >st Y , if E(p{X) > EcpiY) for 
all increasing functions (j) for which the expectations exist. 

Equivalently, X >st Y if and only if P{X > t) > P{Y > t) for all t. 

Definition 2.2 X is larger than Y in the convex sense, X >cx Y, if E(j){X) > E(f){Y) for all 
convex functions (j) for which the expectations exist. 

Note that X >cx Y implies EX = EY and Var{X) > Var{Y). In this sense the convex 
ordering is an ordering of variability in random variables. 

Definition 2.3 X is larger than Y in the increasing concave sense, X >icv Y , if E(j)(X) > 
E(f){Y) for all increasing concave functions (p for which the expectations exist. 

Definition 2.4 X is larger than Y in the Laplace-transform sense, X '>lt Y, if E[e~^^] < 
E[e~^^] for all > for which the expectations exist. 

Note that X >cx Y implies X <icv Y, which in turn implies X <lt Y. 

Finally, for reasons of brevity wc use the following notation. When we say X = [Y \ Z = z], 
we mean that P{X = x) = P(Y = x\Z = z) for all x. 



3 



3 Preemptive LCFS disciplines 



3.1 LCFS preempt-resume 

We first consider the queue with the LCFS preempt-resume (LCFS-p-resume) disci- 

phne. That is, the customer that has been in the system the least amount of time is always 
served, and newly arriving customers preempt earlier arrivals already in service. Within a batch 
customers are arbitrarily labeled, so that we may think of them as arriving sequentially, though 
immediately after each other. Thus, one customer in a newly arriving batch will be considered 
the most recent arrival and will immediately enter service, and the rest of the batch cannot be 
served until that cTistomer, as well as all customers arriving in later batches that preempt that 
customer, are served. 

Customers who resume service after being preempted start their service where they left 
off. Hence, a random service with service time S that is preempted when t units of service 
have already been received has remaining service time [S" — > t]. We also assume service is 
non-idling. Let T be a generic interarrival time, where T has an exponential distribution with 
rate A, and let X be a generic batch size with arbitrary distribution and mean /x. We assume 
that the queue is stable, X^ES < 1. 

Let customer be the last customer in the first batch in the busy period, i.e., the first 
customer to enter service, and let Sq be the service time of customer 0. Let N = N{So) be 
the number of Poisson batch arrival times that occur during the service of customer 0, and let 
N{s) = [N{So)\So = s]. Note that the service wiU be interrupted if N{So) > 0. Let Xq be 
the number of customers in the first batch of the busy period and define M{k,n) = [M\Xq = 
k,N = n] and M{k) = [M\Xo = k], so that M{Xo,N) = M = M{Xq). Let M^, i = 1,2, be 
i.i.d. copies of M, and define maxi=i^...^„ Mj to be if n = 0. For the LCFS-p-resume discipline 
we then have the following characterization of M{k,n). 

Theorem 3.1 The maximum queue length M{k,n) for the /G/1 queue under the LCFS- 
p-resume discipline satisfies 

M(fc,n) =d max{ max Mj /s;M(fc - 1)}, A; > 1, n > 0, (1) 

i=l,...,n 

where M(0) = 0, and Mi, i = 1,2, . . ., and M{k — 1) are independent. 

Proof. We can think of constructing the busy period, conditional on Xq = k, Sq = s and 
N{s) = n, as follows. Denote the arrival epochs, on a clock that only ticks when customer is 
being served, hy < ti <■■■< tn < s. A batch of customers arrives at time ti and starts a 
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new independent busy period (and stops our clock temporarily), except that there are k more 
customers in the queue (the original customers) throughout that busy period. When this first 
sub-busy period is over, at time ti + r say, then customer returns to service and our clock 
resumes ticking. Another batch arrives at time t2 + t, starting a new independent busy period, 
and so on, until the n sub-busy periods have completed, as well as the original service time s. 
Then a new busy period starts with the other k — 1 customers that arrived in the first batch, 
and the maximum queue length during that busy period has the same distribution as M{k — 1). 
Because the arrival process is memory less, this construction is stochastically equivalent to the 
dynamics of a generic /G/1 busy period starting with k customers. □ 

Let P{k,b) = P{M{k) < 6), and P(6) = P{M < b) = EP{Xo,b). So P(0,6) = 1 and 
P(0) = 0. Using the fact that E[P{b — k)'^] is the z-transform, or probability generating 
function of N evaluated at z = P(b — k), we have from Theorem 13. II that for 1 < k < b, 

P{k,b) = E[P{M + k <b)^]P{k-l,b) 
= E[P{b-kf]P{k-l,h) 
= Ee-^(i-^(^-'=))^P(A:-l,6). 

Corollary 3.2 Forb>k>l, 

k 

P(fe,6)=nEe-^(i-^(''-*))^. 
1=1 

If we restrict ourselves to unit batch sizes only, so X = 1 and P{b) = P(l, b), we have the 
following corollary. 

Corollary 3.3 IfX = l, then for b>l, 

P(6) = £e-^(i-^(^-i))^. 

Now we can see how the distribution of S affects M. 

Theorem 3.4 For M^/G/1 queues operating under the LCFS-p-resume discipline, if S' <lt 
S, then M' <st M. In particular, if S' >cx S, then M' <st M. 

Proof. To show that M' <st M, we show that P{k, b) < P'{k, b) (with the obvious definition 
for P') for all k and b by induction on b and k. For each b we have P(0, 6) = 1 = P'(0, b), and 
P[k, 1) = = P'{k, 1) for k>l. Since S' <lt S, 

P(l, 1) = P{S <T) = Ee-^^ < Ee-^^' = P'(l, 1). 
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Suppose P{i,a) < P'{i,a) for a < 6 and all i > 0, so P{a) < P'{a) for all a < b, and suppose 
P{hb) < P'{i,b) for all < i < fc, and consider b and k. From Corollary 13.21 the induction 
hypothesis, and the assumption S' <lt S, it then follows that 

k k 

i=l i=l 
k 

i=l 

This completes the proof. □ 

This theorem is illustrated by Figures Ql^l and ISbelow, showing the probabilities P(M < n) for 
(convexly ordered) families of uniform, Pareto and hyperexponential distributions. 

Remark It is well known (Kelly, 1979) that the M/G/1 queue under the LCFS-p-resume 
discipline exhibits service time insensitivity in the sense that the marginal distribution of the 
number in the stationary system, L, depends on the service-time distribution only through its 
mean. At first this seems at odds with our results, but we must bear in mind that the maximum 
number in the system during a busy period depends on the sample-path evolution of the queue 
length over a busy period. Hence the behaviour of M and L may be very different. This idea 
is further illustrated by the following heuristic example. 

Example Let M be the maximum number in system for an M /G/1 LCFS-p-resume 
queue with 5 = 1 (call this system 1) and let M' be the corresponding maximum when the 
first service time in a busy period, S", is equally likely to be e or 2 — e so S" <cx S', and the 
other service times in the busy period are identically equal to 1 (call this system 2). Then, for 
£ very small, the first busy period in system 2 is equally likely to be very short and have a 
maximum of 1, or it will essentially consist of two busy periods, each evolving as a busy-period 
in system 1. The second of these busy periods starts when the initial customer has received 
1 — e/2 service. That is, roughly, M' is equally likely to be 1 or to have the same distribution as 
max{Afi, M2}, so M' M. Note however that L and L' have roughly the same distribution. 
Indeed, P{L = 0) = P{L' = 0), since the workload is the same in both systems. Furthermore, 
a random arrival during a busy period in system 2 will either see a customer with 5' = e in 
service, with very small probability, or will arrive during one of the two busy periods that each 
evolve as in system 1. Hence L' and L have roughly the same distribution. Finally note that 
the distribution of the length of a busy period does depend on the distribution of S. □ 

The /G/l/b LCFS-p-resume queue also exhibits insensitivity, i.e., the distribution of 
the number in system, Lf,, depends on the distribution of S only through its mean. Hence, the 
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loss rate in the M^/G/1/6 queue, P{Li, = 6), is insensitive to the distribution of S. In contrast, 
our result shows that the probability of at least one loss during a busy period, P{M > 6), does 
depend on the distribution of S, and is greater when S is larger in the Laplace-transform sense. 



3.2 LCFS preempt-repeat with resampling 

Now we suppose that when services are preempted they must be restarted from scratch. The 
new service time is assumed to be an independent random variable with the same distribution. 
We call this the LCFS-p-repeat (with resampling) discipline. Of course, the behavior of the 
queue under the LCFS-p-resume and LCFS-p-repeat disciplines is the same when service times 
are exponential. 

We use the same notation as in the previous subsection. Now, for stability, we need 
XlJ.ESe{S) < 1 and \fiESe{S') < 1, where Se{S) is the effective service time, i.e., the total 
time a random customer must spend in service, including restarts due to interruptions. Thus, 

ESe{S) = E{S A T) + P{S > T)ESeiS), 

where a Ab = min{a, b}, and hence 

ESe{S) - p^g^y (2) 

For T exponential with rate A, it is not hard to show that 



1 - E{e-^''^'^ 
XE{e- 



= ^-F773a5^ (3) 



and hence for stability we need E{e > -|- 1). 

For the LCFS-p-repeat queue, we can identify the following embedded random 

walk. The number in the system at arrival and departure epochs during a busy period is 
equivalent to a random walk on the nonnegative integers with absorbing state 0. The random 
walk starts at the random point Xq, decreases by 1 if T > 5 (a departure), and increases if 
T < S (an arrival). When it increases, it increases by X, where X is independent of S and T. 
Thus, we have the following characterization of M, where I = 1 if T < S" and otherwise, and 
other definitions are as in previous sections. 

Theorem 3.5 The maximum queue length M{k) for the queue under the LCFS-p- 

repeat discipline satisfies 

M{k) =d IM{k + X) + (1 - /) max{k, M{k - 1)}, 
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where M(0) = 0, and I, X, and M{k — 1) are mutually independent, and M{k + X) is inde- 
pendent of I and M{k — 1). 

Let J' be 1 if T > 5", and otherwise. If P(T > S') > P{T > S), then /' >st L Prom 
Theorem 13 . 51 and a couphng argument it then follows that M' <st M. Therefore, we have the 
following. 

Theorem 3.6 ForM^/G/1 queues operating under the LCFS-p-repeat discipline, if E{e ^^') > 
E{e~^^), then M' <st M. 



Note that for the LCFS-p-repeat discipline, we only need for the Laplace transform of the 
service time evaluated at (the arrival rate) A to be ordered for two service-time distributions, 
rather than a complete Laplace-transform ordering. Thus, all possible distributions of service 
times can be completely ordered, and hence we have a complete stochastic ordering of the 
corresponding maximum queue lengths. Of course, it is also true that 5' >ct, S implies S' <lt 
S, which in its turn implies E{e~^^ ) > E{e~^^). 



3.3 LCFS preemptive repeat without resampling 



For our final model, we suppose again that when services are preempted they must be restarted 
from scratch, but now the service time is only drawn from the service-time distribution once. 
We call this the LCFS-p-repeat-without-resampling discipline. Note that the LCFS-p-repeat 
and LCFS-p-repeat-without-resampling disciplines are the same for deterministic service times. 
For stability, we need again XfiESe{S) < 1 and XfiESe{S') < 1, where Se{S) is the effective 
service time. Given S = s, the service time is deterministic and the effective service time Se is 
the same as in equation (jj}, that is 



Hence, 



ESe = ESsiS) = E{E[Se{S) \S]) = ^[^e^^ - 1]. 

A 



So for stability we need E{e^^) < + 1)/^. If, for example, S is exponentially distributed 
with mean z^, then for stability we need v > X{fi + 1). Note that this value is larger than for 
the repeat-with-resampling discipline. Intuitively, a large value of the service time has a large 
probability of being interrupted and having to start over, and each time it restarts it will again 
have a large service time. 
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With Xq, Sq, and M defined as in the last subsection, and with Ti defined to be the first 
interarrival time after the busy period starts, we now let M{k,s) = [M\Xq = k,SQ = s] and 
M{k) = M{k, S) = [M\Xo = k]. Let I(s) = 1 if Ti < s and otherwise. We have the following. 



Theorem 3.7 The maximum queue length M{k,s) for the M^/G/1 queue under the LCFS- 
p-repeat-without-resampling discipline satisfies 

M{k, s) =d I{s) max{M + k; M{k, s)} + (1 - I{s)) max{k, M{k - 1)}, 

where M(0) = 0, and where I{s), M, M{k,s), and M{k — 1) are independent. 

Proof. Given Xq = k and 5o = s, if an arrival occurs before the first service completion a 
new i.i.d. (sub-)busy period starts, except that there are k additional customers in the queue. 
When that sub-busy period ends, the original busy periods starts again, independently of Ti and 
of M for the ending sub-busy period, with Xq = k and 5o = s. If the first service completes 
before an arrival, then we may consider the remainder of the busy period as a new, independent 
busy period with k — 1 initial customers. □ 

Let P{k,b,s) = P{M{k,s) < b), P{k,b) = P{M{k) < b) = EP{k,b,So), and P{b) = 
P(M < b) = EP{Xq, b, Sq). We have the following corollary to Theorem 13.71 

Corollary 3.8 For b > k > I and for all s, 

P{T > ,s)P{k-l,b) _ e-^'Pik -l,b) 



P{k,b,s) 



1 - P{T < s)P{b - k) l-il-e-^')P{b-k)' 
Using this corollary, we can show the following. 

Theorem 3.9 For /G/l queues operating under the LCFS-p-repeat-without-resampling 
discipline, if S' <icv S, then M' <st M. Hence, if S' >cx S, then M' <st M. 

Proof. From our corollary above, for b > k > I, and s > 0, 

P{k-l,b) 



P{k,b,s) 



e^'il - P{b - k)) + P{b - k) 



It is easy to show that f{s) := a/{ce^^ + d) is a decreasing convex function of s for all a,c,d> 0, 
c + d > (so — /(s) is increasing and concave). Hence, if S' <icv S then —Ef{S') < —Ef(S) 
and Ef{S') > Ef{S). Also note that P{k,b,s) is increasing in P{k — 1,6) and P{b — k) for 
fixed s. The result now follows using an induction argument similar to the one in the proof of 
Theorem 13.41 □ 
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4 Numerical Illustrations 



For the LCFS preempt-resume discipline, we calculated P{M < n) using Corollay I3..SI The 
figures below show the results for several convexly ordered families of distributions, illustrating 
Theorem 13.41 



uniform(0.9,1.1) 
uniform(0,2) 



Figure 1: P(M < n),n = 1,2, ...,20, for uniform service-time distributions with ES = 1; 
A = 0.9. 



Pareto{1 .5) 
Pareto(2) 
Pareto(IO) 



Figure 2: P{M < n), n = 1,2,..., 20, for Pareto(a) service-time distributions with distribution 
function Fa{x) = 1 - {{a - > (a - so ES = 1] X = 0.95. 
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hyperexponential(1 ) + 
hyperexponential(3) x 
hyperexponential(1 0) >i 

X =•< X "5 
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Figure 3: P{M < n), n = 1, 2, . . . , 10, for hyperexponential(A;) service-time distributions with 
P{S = Xi) = 1 — = 1 — P{S = X2), where Xi and X2 are exponentially distributed with 
mean 1/(2(1 - 2"^)) and 2^/2 respectively, so ^^S' = 1; A = 0.95. 
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